Analytical expressions of the dispersive contributions of the nuclear optical model 
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Analytical solutions of dispersion relations in the nuclear optical model have been found for both 
imaginary volume and surface potentials. A standard Brown-Rho shape has been assumed for 
the volume imaginary term and a Brown-Rho shape multiplied by a decreasing exponential for the 
surface contribution. The analytical solutions are valid for any even value of the exponent appearing 
in these functional forms. 
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A significant contribution to the optical model theory 
during the last two decades can be considered the work 
of Mahaux and co-workers on dispersive optical model 
analysis Jl], [|, |[ §, ||. The unified description of nuclear 
mean field in dispersive optical model is accomplished 
by using a dispersion relation, which links the real and 
absorptive terms of the optical model potential (OMP). 
The additional constraint imposed by dispersion relation 
helps to reduce the ambiguities in deriving phenomcno- 
logical OMP parameters from the experimental data. 
In a dispersion relation treatment, the real central po- 
tential strength consists of a term which varies slowly 
with energy E, the so called Hartree-Fock (HF) term, 
V HF (r, E), plus a correction term, AV(r,E), which is 
calculated using a dispersion relation. The depth of the 
dispersive term of the potential AV{r, E) can be written 



1 ' W-oo E'-E 



(1) 



Being W(r, E) the imaginary part of the OMP. Assuming 
that W(r, E = E F ) = and AV(r, E = E F ) = 0, where 
E F is the Fermi energy, then equation (0) can be written 
in the substractcd form: 



AV(r,E) = - / W(r,E') 
n J-oo 



E'-E 



Ej- 



dE' 
(2) 

With the assumption that W(r, E) be symmetric respect 
to the Fermi energy, the equation (^|) can be expressed 
in a form which is stable under numerical treatment 
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namely: 

AV(r,E) = -(E-E F 



' W(r,E') -W(r,E) 
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-dE' 
(3) 

The dispersive term AV(r, E) is divided into two terms 
AV V (r, E) and AV S (r, E) , which arise through disper- 
sion relations (|l|) from the volume W v (r, E) and surface 
W s (r, E) imaginary potentials respectively. If the imag- 
inary potential geometry is energy dependent, the radial 
dependence of the dispersive correction can not be ex- 
pressed using a Woods-Saxon form factor f ws (r). How- 
ever, to simplify the problem, the OMP geometry pa- 
rameters are usually assumed to be energy independent. 
In this case, AV(r,E) = AV(E)f ws (r) and the energy 
dependence of the real volume V v (E) and surface V s (E) 
parts of the dispersive OMP are given by: 



V V (E) = V HF (E) + AV V (E) 
V S (E) = AV S (E) 



(4) 



It is useful to represent the variation of surface W s (E) 
and volume W v (E) absorption potential depths with en- 
ergy in functional forms suitable for the dispersive optical 
model analysis. An energy dependence for the imaginary 
volume term has been suggested in studies of nuclear 
matter theory M: 



W V (E) = A, 



(E — E F 



(E-E p )n + (B v y 



(5) 



where A v and B v are undetermined constants. Brown 
and Rho 0, propose n = 2, while Mahaux and Sartor Q 
have suggested n = 4. 

An energy dependence for the imaginary-surface term has 
been suggested by Delaroche et al 11 to be: 



W S (E)=A S 



(e-e f y 



(E-E F )^ + (B S 



exp(-C s |£ - E B 



(6) 
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where m = 2,4 and A S ,B S and C s are undetermined 
constants. 

According to equations (^|) and (^) the imaginary part of 
the OMP is assumed to be zero at E = E F and nonzero 
everywhere else. A more realistic parameterization of 
W v (E) and W s (E) forces these terms to be zero in some 
region around the Fermi energy. A physically reasonable 
energy for defining such a region is the average energy of 
the single-particle states E p B. 



Therefore a new definition for imaginary part of the OMP 
can be written as: 



Wy(E) = 



(E-E F 



(E-E r 



-(B v 



E v 



and likewise for surface absorption. 
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W S (E) = 



A s e 



-C S \E-E p \ 



E. 
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(E-E p )™ + (B S 



< E < E p 
E > E r _ 



(8) 



The symmetry condition 



W(2E F — E) = W(E) 



(9) 



is used to define imaginary part of the OMP for energies 
below the Fermi energy. Equations are used to 

describe imaginary absorptive potential in this contribu- 
tion. 

Nonlocality effects on the imaginary part can be treated 
according to the work of Mahaux and Sartor jl| and the 
analytical solutions of their contributions coming from 
the dispersion relations have been found by VanderKam 
et al. j8|]. 

In a recent work we have presented a general numer- 
ical solution of the dispersion integral relations between 
the real and the imaginary parts of the nuclear optical 
potential Q. In this work we develop analytical so- 
lutions of dispersion relations, equation (^|) for usually 
employed functional forms of the imaginary potential, 
W(E), used in dispersive optical model analyses (equa- 
tions (f7}|j|)). According to Q we introduce the quantities 
E = E p — E F (the offset energy), E + = E x + E and 
E- = E x — E . U and E x are the integration variable 
and the excitation energy (E — E F ), respectively. In the 
next formulae, we will drop the subscripts on the con- 
stants A, B, C. 

For the surface potential W S (E), given by equations (||) 
and (||), we can write the dispersive integral (^) for m 
even as 



E f°° 
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E)(E'-E F 
U m exp(-CU) 
u (U m + B m )(U - E-)(U + E ) 
U m cxp{~CU) 
(U m + B m )(U + E+)(U + E a ) 



-dE' 



dU 



dU (10) 



The integrand can be replaced by its expressions in terms 
of poles and residues: 



E. r 



U" 



tt (U m + B m ){U T E T )(U + E ) 



1 }™ Re S ( Pj ) Res( T E ± ) Res(E ) 



where the pj are the in zeroes of (U m + B m ) and the 
Res{pj) their residue, that is : 



Pj = B exp I i ■ 7r 



Res{pj) 



Er 



m 
Pj 



m ( Pj T E T )( Pj + E ) 



(12) 
(13) 



where and — E are the poles of U =F E T and U + 

E , whereas Res(±E^) and Res(—E ) are their residues 
respectively 



Res(±E T ) = ± 



Res(-E ) = =f 



(E T ) m + B r - 

(E ) m 
{E ) m + B r ' 



(14) 
(15) 



As was pointed out by Raynal |l0| the contribution of 
each complex pole pj to the surface dispersive integral is 



Res(p 3 )e- cu 



U-pi 



Res(pj)e 



-c Pj 



dU 



ex P( — z ) 



dz 



= Res(pj)e- Gp * Exi-Cpj) 



(16) 



and the same expression without exponential for the vol- 
ume potential W V {E), given by equations (^) and (^). 



where Ei(z) is the Exponential Integral Function E\ pi) . 
For the real poles corresponding to the second term in the 
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right side of (p_lj), the contribution of each one is given by: 



Res(E-) e- CE -Ei(CE_) 



(18) 



Res(TE ± )e- cu 



Res{TE ± )e TCE TV 



-dU 



Res(TE ± )e TCE ^ Ei(±CE ) 



(17) 



where Ei(x) is the Exponential Integral Function Ei [ p"T| . 
The contributions coming from the third term of the rigth 
side of eq. ([ll]) in the integral (|Io|) cancel. 
For the volume dispersive integral the exponential is 
missed in the integrand, therefore each pole gives a di- 
vergent contribution. Nevertheless their sum is a finite 
quantity, which can be calculated by taking the proper 
limit. 

We are quoting below exact expressions for the surface 
and volume dispersive integrals for any even values of m 
and n in the surface and volume imaginary potentials re- 
spectively. This make a difference to ref. M, where, in 
surface dispersive integral with m = 4, approximations 
were made limiting the energy range of validity for given 
analytical expressions. 

For the case of the dispersion relation using surface imag- 
inary potential W s (E) with m even according to equa- 
tions (|) and (0), the dispersive correction is 



where Zj comes from the sum of residues Res(pj) in the 
two integrals ( |l0| ) and is given by 



= E x P] (2 P] + E + - E_) 
3 m ( Pj + E )(p s + E + )( Pj - E_) 



(19) 



For the case of the dispersion relation using the volume 
imaginary potential W v (E) with n even according to 
equations (0) and ([|), the dispersive correction is 

A f " 

AV V (E) = -- |X>In(-Pi) 
+ Res{-E + ) InE^ 



+ Res(E-) ln\E. 



(20) 



In conclusion, we have developed analytical solutions 
of dispersion relations for the volume and surface terms 
of the OMP. The formulae are compact and easy to 
implement in current generation of the optical model 
parameter search codes. We stress that the solutions are 
valid for any even exponent norm. 
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AV S (E) = -^Z j e-^ c E 1 {-p j C) 
- Res(-E+) e CE+ Ei(-CE+) 
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